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Materials and Methods 


Finite-element calculations for simulating non-ideal materials 


Because of the employed superconductor and ferromagnetic materials, some non-ideal behavior would 
be expected. In order to quantify these effects we have performed finite-element calculations under the 
following assumptions: 

- The structure is infinite along its cylindrical axis (this means two-dimensional distributions of cur- 
rents and fields have been computed). 

- Superconductor behaves according to the critical state model (27) - this approximation is widely 
used for high-temperature superconductors with strong pinning of magnetic flux working at 77 K in 
quasi-static conditions (frequencies below | kHz). 

- Magnetic permeability of ferromagnetic material is finite and isotropic. 

The calculation procedure has been developed earlier for the purpose of investigating the possibility 
to improve the transport capability of a superconducting wire with the help of a ferromagnetic layer 
(28). It provided excellent predictions for the electromagnetic properties of high-temperature tapes with 
ferromagnetic substrate (29). As an example Fig. S1 shows the distribution of magnetic field in the 
magnetic cloak shell with number of layers nsc = 2 and npm =7 at switching on the uniform magnetic 
field of 40 mT. Magnetic invisibility of the magnetic cloak is evident - a deformation of straight induction 
lines is hardly observable. 

Quantitatively, the effect of combining SC and FM materials is illustrated by the calculated results in 
Fig. S2 showing the scan of the vertical component of magnetic induction By, at a fixed height above the 
two-shell cloak (along the dashed line in Fig. S1). 

The properties of materials entering the calculations have been derived from independent tests. From 
the critical current of the 12 mm wide tape measured at 77 K the critical current density in a super- 
conducting layer of supposed thickness 1 jum was evaluated as Je = 3 x 10!° A/m?. This property 


characterizes the transport capability along the tape length i.e. in circumferential direction in our geom- 


etry shown in Fig. S1. When putting this value in our simulations the results predicted much stronger 
screening than observed experimentally. Eventually, nice agreement was achieved when assuming Je = 
10!°A/m?. Such difference can be explained by the combination of several factors. First, the currents 
securing the cloaking effect flow in direction perpendicular to the drawing plane of Fig. S1 (along z- 
coordinate) that is across the tape width and commercial tapes are not optimized for this direction of 
current. Even more significant is the fact that all the screening current must change direction at the ends 
of a finite tubular cloak. This flow of current - particularly limiting in our case of L comparable to Qin - 
happens along the tape edges, well known for being the place of a reduced critical current. We believe 
these end effects cause the observed reduction of transport capability by factor 3 compared to a straight 
tape. 

The magnetic permeability of the FeNiCr tape was measured on a toroidal sample made of 3 turns. 
In wide scale of excitation fields the relative permeability ranged between 11 and 18. The difference in 


predicted distributions obtained assuming these limiting values was negligible. 


Construction of the bilayer 


A 1 um thick layer of cuprate superconductor ReBCO deposited on 100 wm metallic substrate (20) is 
able to carry without resistance the electrical current exceeding 400 A in 12 mm wide tape when cooled 
by liquid nitrogen to the temperature of 77 K. Changing the number of turns the tape has been wrapped 
around the central bore, nsc, allows to vary the strength of magnetic field expulsion. An outer shell with 
the task of sucking the magnetic flux expelled from central bore is wound from 100 um thick Fe1gCrgNi 
alloy sheet. With the number of FM turns, npm, we can adjust the magnetic response of the outer shell. 
Thickness of this tape is another free parameter of the geometry. 

On the cylindrical mandrel from textile-reinforced plastic two turns of superconducting tape 12 mm 
wide were first wrapped. To avoid uncontrolled short-circuiting of the turns, 70 um thick kapton layer 


was inserted between them. Secondly, from the sheath of FeNiCr alloy a strip 12 mm wide was cut in 


the length sufficient to wrap seven turns of this ferromagnetic material, alternated with a 215 um thick 
kapton layer placed between them. 

This exterior ferromagnetic part, consisting in an alternation of ferromagnetic layers (with measured 
permeability between 11 and 18) and kapton layers (with null magnetic contribution), would yield a 
decreased effective permeability much closer to the theoretical value of 2 = 3.5 calculated with Eq. 1 


using the final radii relation R2/Rı = 1.34 of the experimental set-up. 


Preparation of the experimental bilayer 


Both the high temperature superconductor and the ferromagnetic alloy used to manufacture the two shell 
structure exhibit some magnetic hysteresis. This means the magnetic profiles measured after various 
histories of applied field could differ. In the case of superconductor the remanent currents could be 
canceled by warming-up above the critical temperature. For the ferromagnetic material one can achieve 
a quasi-demagnetized state after application of a series of magnetic field pulses with alternating polarity 
and decreasing amplitude. We used both these measures to obtain in each of the above mentioned cases 
the profile corresponding to the history of applied field that is a simple change from zero to the desired 


value. 


SOM text 


Analytical derivation of cloaking behavior for a superconductor-magnetic cylindrical bilayer 
in uniform applied field 


We consider a cylindrical shell (infinite in the z direction) of radii Rọ and Ry > Ro with uniform 
permeability jz; and another concentric cylindrical shell of radii Ry and Rə > Rı with permeability u2. 
A uniform field H, is externally applied in the y direction. 

We want to calculate the magnetic field distribution in the whole space. Since no free currents are 


present, the magnetic field H can be written in terms of a magnetic scalar potential ¢, H = —V 4, that 


should satisfy the Laplace equation V7¢ = 0 in all the regions of space. For the present geometry, 
the general solution of Laplace equation expressed in cylindrical coordinates, where @ is the angle with 
respect to the applied field direction, is (30) 


CO 
$i = aoi + boilnp + 5 an ip” sin(n + ani) 
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a 5 bnip ” sin(nO T Bni), (S1) 
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where an, bn, Qn, Bn are constants to be determined by the boundary conditions and the subindex i stands 
for the different regions where ¢ is going to be determined: 7 = 0 for the region inside the hole (p < Ro), 
i = 1 for the inner layer (Ro < p < Rı), i = 2 for the outer layer (Ry < p < Rə), andi = 3 for the 
exterior region (p > Rə). 

Taking into account that (a) ¢3 should tend to — Hap cos 0 (the applied scalar potential) when p — oo, 
(b) ġo should be finite at p = 0, and (c) the standard magnetostatic boundary conditions must be obeyed, 


then one can find the scalar potential in all the regions. The final solutions can be expressed in the form 


ġo = Fpcosé, (S2) 

og = (Dp + =| cos 0, (S3) 

g2 = (Bo + <) cos ð, (S4) 
A 

o3 = (tp + 2) cos 0, (S5) 


where A, B,C, D, E, and F are the solutions of the following system of equations 


BR?+C-DRÌ-E = 0, (S6) 

DR+ E-FR? = 0, (S7) 
-A-umBR3+mC = HR}, (S8) 
A-BR2-C = H,R3, (S9) 
2BRi—j2C-j4DRi+mE = 0, (S10) 


mDR — mE- FR = 0, (S11) 


which can be analitically obtained. The general combination of positive permeabilities for the cylindrical 


case that yields a zero distortion of the applied field in the region p > Ro is 
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where we have defined an auxiliary parameter 

K = (uj — 1)(Ro — RI)(RĪ + R3). (S13) 


AS a particular case, when the inner layer is superconducting (11 = 0), the previous constants reduce 


(1+ po) R7 — (—1+ po) R5 
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F = 0. (S17) 
In this case, by setting A = 0 in Eq. (S14) we obtain Eq. 1 of the main text 
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Superconductor-magnetic spherical bilayer in uniform applied field 


Analytical derivations can also be made for a spheric configuration, instead of a cylindrical one. A hollow 


sphere with a hole of radius Rọ, an inner layer with radii Rọ and Rı and permeability u1, and an outer 


layer with radii Rı and Ro with permeability u2 (Ro < Rı < Rə) produces no external field distortion 
when 


P 
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where we have defined another auxiliary parameter 


P = (1 +443) RY + (2+ p) (1 +241) RE RS — (—1 + wi) RG((1 + 4m) R] +2(1 + 1) R3). (S20) 


Magnetic fields in all regions can be obtained analytically also in this case. When the inner layer is 


an ideal superconductor, the non-distortion in the exterior regions requires 


2R3 + RÌ 
SS $21 


Although not directly related to our work, some derivations for the electric field in systems with 
coated spherical inclusions were presented in the context of the theory of composites (31,32). There, the 
goal was not to design cloaks but to find the effective conductivity for a composite comprised typically 
of spherical grains that did not distort an applied electric field. In (33), a discussion of the required 
dielectric constant values for a compound dielectric ellipsoids to reduce its scattering of electromagnetic 


radiation was presented. 


Superconductor-magnetic cylindrical bilayer in non-uniform applied field 


Here we discuss the cloaking behaviour of a superconductor-magnetic cylindrical bilayer when the ap- 
plied field has some spatial variation. We provide a way to estimate the external field distortion and a 
general procedure to reduce it. The distortion zones showing where the magnetic induction B differs 
from the externally applied magnetic induction Bext in more than p per cent are calculated as the points 


in space where 
| B i Bext | 


> p/100. S22 
Ba] p/ (S22) 


The magnetic response of a superconducting-magnetic bilayer to the fields of a single small mag- 
net (basically a dipole field) at different distances from the bilayer are shown in Fig. S3, calculated by 
Comsol Multiphysics software, using the electromagnetics module (magnetostatics). In these cases of 
non-uniform applied field, a distortion of the external field occurs - distortion gets larger with increasing 
applied field inhomogeneity-, but it can be seen that this distortion is less than few percent in all space 
except in a region close to the field source and the bilayer. As an example, the dark pink region corre- 
sponding to a 5% distortion that in the case of a dipole at a distance 2.52 occupies a region of about 
2-3 Ro, extends to more than 15 Rə in the case of a single magnetic layer or a single superconducting one 
(like those in Figs. 1, A and B). 

There is an important effect worth to discuss. When reducing the thickness of the outer magnetic 
layer Ro, keeping Rı constant and changing u2 according to Eq. 1, the region in which the external 
field is distorted is reduced. Actually, for thin magnetic layers, Rə — Ri << Rı, Eq. 1 reduces to 
u2 ~ Rı/(Rə— R}). In that case, the large values of u2 allow a strong modification of the external field 
inside the thin layer. This effect is seen in Fig. S4, where the calculated distortion zone in the case of 
a near dipole, an unfavorable case in terms of field inhomogeneity, is clearly reduced by decreasing the 


thickness of the magnetic layer. 


Figure S1: Example of 2D calculation for uniform magnetic field applied in perpendicular direction to 
the two-shell magnetic cloak (that is assumed infinite in the direction perpendicular to the drawing plane). 


Two inner turns are from superconducting material (thickness 1 zm), seven outer turns (thickness 100 
um) from ferromagnetic alloy. 


B (mT) 
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—— FM only 
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25 1 L 1 li . 1 1 
-10 -5 0 5 10 
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Figure S2: Calculated distribution of y-component of magnetic induction along the dashed line shown 
in Fig. S1 for applied field 40 mT, superconductor critical current density 10!° A/m?, permeability of 
magnetic material up = 11 (blue line). For comparison the calculated results are shown for the case of 
switched off superconductor (FM only) and non-magnetic ferromagnetic layers (SC only). 
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Figure S3: Calculated field lines for the interaction of a bilayer device (with Ro/R, = 1.25) with the 
field created by a magnetic dipole placed at different distances from the center of the device: (A) 2A), 


(B) 2.51, (C) 3Rı and (D) 3.51. Bright (dark) pink-shaded regions indicate where distortion is larger 
that 1% (5%). 
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Figure S4: Calculated field lines and distortion regions for the interaction of a dipolar field (source is at 
a distance of 2.5, from the center of the bilayer) with bilayer devices with different thickness for the 
outer layer [(A) R2/Rı = 1.15, and (B) R2/Rı = 1.025], with corresponding permeabilities given by 
Eq. 1. Bright (dark) pink-shaded regions indicate where distortion is larger than 1% (5%). 
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